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ABSTRACT 


FOREWORD 


The  computer  program  described  in  this  report  was  prepared  in  the 
Ballistic  Sciences  Branch  of  the  Computer  Programming  Division,  Compu¬ 
tation  and  Analysis  Laboratory,  by  Miss  Beverly  A.  Cooper.  The  work 
was  performed  under  the  authority  of  NWL  Technical  Assistance  Request 
for  Project  No.  T-06AD  by  the  Computation  and  Analysis  Laboratory  in 
support  of  the  Warhead  Supporting  Research  Program  being  conducted  by 
the  Warhead  and  Terminal  Ballistics  Laboratory. 
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INTRODUCTION 


The  U.  S.  Naval  Weapons  Laboratory  is  currently  preparing  a  hand¬ 
book  series,  "Warhead  and  Terminal  Ballistics  Handbook,"  in  conjunction 
with  the  Warhead  Supporting  Research  Program  established  in  accordance 
with  Bureau  of  Naval  Weapons  Weptask  No.  RMMO-42-G03/210-1/F008-08-06. 
The  objective  of  this  handbook  is  to  provide  systematic  documentation 
of  kill  mechanism  technology,  target  vulnerability  data,  and  warhead 
terminal  ballistic  performance  data.  To  be  included,  is  a  section  on 
the  terminal  ballistic  performance  characteristics  of  bombs,  projec¬ 
tiles,  and  rocket  and  missile  warheads  of  current  tactical  interest  to 
the  Navy.  A  preliminary  edition  of  this  section,  published  as  NAVWEPS 
Report  No.  7673,  is  currently  being  revised  to  include  data  recently 
acquired  through  an  extensive  warhead  test  program  conducted  at  the 
U.  So  Naval  Weapons  Laboratory. 

A  significant  portion  of  this  program  was  a  series  of  arena  tests 
of  statically  fired  warheads  to  establish  measurements  of  the  fragmen¬ 
tation  characteristics --velocity,  density,  mass,  and  spatial  distribu¬ 
tion  of  fragments — for  each  of  several  warheads.  Typically,  these 
tests  are  performed  in  a  field  arena  similar  to  the  one  described  in 
Figure  1.  The  arena  consists  of  two  180-degree  circular  sectors.  One 
sector  is  constructed  from  22-gauge  mild  steel  plate,  and  the  other  is 
constructed  of  celotex  panels.  The  warhead  to.be  tested  is  positioned 
on  a  stand  at  the  center  of  the  arena  with  its  axis  horizontal  and 
intersecting  the  sides  of  the  arena  where  the  two  sectors  meet,  with 
the  forward  end  of  the  warhead  aimed  at  the  O-degree  position  on  the 
wall.  A  ricochet  fence  is  built  about  the  stand  between  the  warhead 
and  the  arena  to  prevent  fragments  that  strike  the  ground  from  reaching 
the  wall  panels . 

i 

High-speed  motion  picture  cameras  are  placed  about  the  steel  sec¬ 
tor  to  record  impacts  of  fragments  from  the  detonated  warhead  on  the 
arena  plates.  A  kilocycle/second  timing  trace  is  superimposed  on  the 
film  in  each  camera  during  operation  to  allow  for  the  determination  of 
the  time  lapse  from  detonation  to  impact.  Given  the  resulting  time 
data  and  the  distance  from  the  warhead  to  the  arena  plate,  an  estimate 
of  the  average  velocity  of  the  fragments  can  be  obtained.  Fragments 
impacting  on  the  celotex  sector  are  recovered  upon  completion  of  the 
test,  weighed,  and  categorized  by  mass  interval  to  obtain  estimates 
of  the  fragment  number  and  mass  distribution  associated  with  the  war¬ 
head. 


A  spatial  distribution  of  the  fragments  is  constructed  by  tabu¬ 
lating  the  above  data  for  each  five-degree  polar  zone  (angular  inter¬ 
val  measured  from  horizontal  axis  of  warhead) . 
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O’  CAMERA 


FIGURE  I.  TYPICAL  ARENA  SET-UP 


DATA  REDUCTION 


To  facilitate  and  expedite  reduction  of  the  raw  test  data  to  a 
useful  and  meaningful  form,  a  high-speed  computer  program  was  formu¬ 
lated  and  coded  in  FORTRAN  IV  for  execution  oh  an  IBM  7030  computer. 

The  flow  chart  given  as  Figure  2  outlines  the  major  features  of 
the  program.  For  each  mass  interval  (weight  group)  within  a  polar 
zone,  the  fragment  density  and  average  mass  are  computed.  Upon  com¬ 
pletion  of  these  computations,  initial  fragment  velocity,  total  frag¬ 
ment  weight,  and  number  of  fragments  per  steradian  for  the  polar  .zone 
are  calculated.  Each  zone  is  considered  in  succession  until  the  total 
fragment  space  has  been  spanned.  These  output  are  then  tabulated  as 
shown  in  Figure  A-3. 

a.  The  fragment  density  for  each  mass  interval  is  obtained  by 
scaling  the  average  number  of  fragments  recovered  by  the  value  com¬ 
puted  from  the  equation: 


fk 

ak 


21T  (cos  -  cos  8^) 


where , 

H 

arc tan  x  =  ^ 

r  =  radius  of  mass  arena,  feet 

H  =  height  of  mass  arena,  feet 

8^  =  angular  distance  from  horizontal  axis  of 
warhead  to  the  lower  bound  of  k-th  polar 
zone 

The  average  mass  is  simply  the  total,  weight  of  the  recovered  fragments 
divided  by  the  number  of  fragments. 
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ENTER 


Figure  2.  Flow  Chart 


ttm&tmxs.'tirri*  wra 


b.  The  initial  fragment  velocity  for  the  polar  zone  is  determined 
from  the  equation, 


V 


V  = 

V  = 
e  = 

and, 

w  = 

for, 

/a  = 

Cd  = 

R  = 

The  value  of 


where, 


b  = 


c 

m 


initial  fragment  velocity,  feet/second 
average  fragment  velocity,  feet/second 
base  of  natural  logarithm 

0.489  /9a<W4} 


o 

air  density,  pounds/ feet 

drag  coefficient 

radius  of  velocity  arena,  feet 

area  mass  ratio,  centimeters^/grams . 

the  area  mass  ratio  is  given  by  the  relationship, 


fragment  area  mass  constant 
fragment  area  mass  exponent 

fragment  mass,  which  in  this  case  is  the  average 
fragment  mass  taken  over  the  entire  polar  zone  . 
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c.  The  total  fragment  weight  per  polar  zone  is  obtained  by  merely 
summing  the  product  of  the  number  of  fragments  and  the  average  fragment 
weight  for  each  mass  interval  over  the  total  polar  zone; 

d.  The  steradian  measure  of  the  polar  zone  is, 


=  27T(cos  0fc+i  -  cos  0^) 


where  0^  is  as  defined  in  subparagraph  a,  above.  Dividing  the  total 
number  of  fragments  in  the  polar  zone  by  Afc  yields  the  number  of  frag¬ 
ments  per  steradian  for  the  zone. 

The  remainder  of  the  pertinent  fragmentation  data — heaviest  frag¬ 
ment  recovered  and  the  number  of  hits  (fragment  impacts  on  velocity 
arena  wall)  from  which  the  average  velocity  is  calculated — is  obtained 
by  observing  the  raw  test  data.  These  values  are  inserted  manually  in 
the  spaces  provided  in  the  output  format.  The  total  process  is  then 
repeated  for  each  successive  polar  zone  until  the  total  fragment  space 
has  been  spanned. 

A  listing  of  the  program  deck,  the  input  format,  and  a  sample  of 
the  input  and  output  from  the  program  are  included  in  Appendix  A. 


8 


APPENDIX  A 


PROGRAM  DECK 


A  listing  of  the  FORTRAN  IV  program  deck  designed  for  execution 
on  the  IBM  7030  is  given  as  Figure  A-l.  / 

INPUT  FORMAT 


A  description  of  the  input  format  is  shown  in  Figure  A-2,  and  the 
various  input  data  required  by  the  program  identified  below: 


Card  Type  1 

NRUNS 

Number  of  fragment  data 
tables  to  be  prepared 

Card  Type  2 

INT 

Number  of  mass  intervals 
per  polar  zone  (■£  20) 

NZONE 

Total  number  of  polar  zones  ( 636) 

MULT 

MULT  =  436/NZONES  (Integer) 

Card  Type  3 

TAB 

Table  number 

WEAP 

Weapon  identification 
(center  in  these  columns) 

Card  Type  4 

EXPLO 

Type  of  explosive 

REFNO 

Reference  number 

NR 

Number  of  rounds  tested 

RADM 

Radius  of  mass  arena 

R 

Radius  of  velocity  arena 

H 

Height  of  arena 

Card  Type  5 

RHOA 

Air  density 

CD 

Drag  coefficient 

B 

Fragment  area/mass  constant 

C 

Fragment  area/ mass  exponent 

Card  Type  6 

NFRAG 

Number  of  fragments  per  mass 
interval  per  zone 

Starting  with  the  first  polar  zone,  list  the  number  of  fragments  in 
the  largest  mass  interval,  the  next  largest  mass  interval,  etc.,  until 
all  mass  intervals  have  been  accounted  for.  Continue  listing  in  the 
same  order  for  each  successive  polar  zone  until  all  polar  zones  have 
been  considered. 
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PROGRAM  TO  EDIT  FRAGMENTATION  DATA 

D I  KENS  I  ON  NFRAG  ( 20  *36  J  »  W(  20*36  )  » V3AR  t  36. )  *  THETA  ( 37 ) , AK( 37 ) »  CAK  ( 3-7  ) 
1 * CNBAR  ( 2Q*,6  )  »C.KBAR  ( 20*36 )  *NO(36)  ,  VOI-36  )  ,  AINTI  20  )  ,RAp(  37 )  » 

2WEAPI 1C)  *EXPLO(  5  )  »  FRAG! 20*36 )  *REFNC(  2  )  ,ASIN(37)  ,WBAR(36) 
EQUIVALENCE  (NFRAG.FRAG) 

10  FOR'-'AT  ( 1 1 A5  ) 

’1  F0PVAT(7A3,I5)3F5.1 ) 

12  FORMAT (4F]C. 5) 

13  FORMAT (141 5) 

14  FORMAT (  7F10.5) 

20  FORMATt  1H1 , 44X * 6HTABLE  »  A5 » 18HFRAGMFNTAT I  ON  DATA/37X * 10A5 ) 

21  FORMAT! ] 1H1AVERA6E  0F,I5»/H  ROUNDS /22H  RADIUS  OF  MASS  ARENA  ,F5.1* 
14H  FT./26H  RADIUS  OF  VELOCITY  ARENA  »F5. 1 »4H  FT./19H  HEIGHT  OF  CEL 
2CTEX  ,F4.1*4H  FT./16H  TYPE  EXPLOSIVE  ,5A5/15H  REFERENCE  NO.  ,2A5> 


22  FORMAT  ( 1 26H.K  ,  - - - - - - - 

1----  -  -  -  -  FRAGMENT  WEIGHT  GROUPS  (GRANS)  -  -  -  -  - 

2  -  -  ) 

23  FORMAT  ( 13H  POLAR  *  F5.1,2H  •»  ,3X»6!4XF6.3  *  1H-,  F6. 3  )  ) 

24  FORMAT ( 123H  ZONE  N  M  N  V  N  v 

1  N  M  N  v  N  (•’  N 

2  M) 


25  FORMAT!  1H  *  13  * 1H-* 1 3 *2 ( 1XF7.3 ) ,6 < 2XF7.3 , 1XF7. 3 ) ) 

26  FORMAT (52HON=AVERAGE  NUMBER  OF  FRAGMENTS  PER  TOTAL  dolaR  ZONE./ 
133H  M=AVER AGE  FRAGMENT  MASS  (GRAMS). /46H  V=D1 STANCE ’ OF  FRAGMENT  Fl 


2  I ght /TIME  OF  FLIGHT.) 

27  FORMAT (1H5,//) 

28  FORMAT  (50HL - ,38X, 

122HTOTAL  NO.  .OF  FRAG.)' 

29  FORMAT!  2H  *2 ( F6. 3 » 1H- *F6 .3 *4X ) >F6. 3  * 1H-»F6. 3 »56H  HEAVIEST  V 

1EL0CITY  (FT/SEC)  FRAGMENT  PER/108H  N  M 


2N  M  MM  FRAGMENT  HITS  AVERAGE  INITIAL  WE 

3  I GHT  ( GMS )  STERADIAN) 

30  FORMAT  ( 1H  , F7 . 3 » 1 X » F7 . 3 » 2 < 2X , F7. 3 ♦ IX * E7 . 3 ) »  16X. 2 ( 1 XF7.2 ) » 

13XF9. 2.7X15) 

3 1 ' FORMAT ( 1H.  >84XF9,2) 

READ  13.NRUNS 

THETA( 1 ) =0.0 

RAD< 1 ) =0.0 

ASIN(1)=0.0 

DO  400  LL=1 *NRUNS 

READ  13* INT.NZONE* IMULT 

READ  10.TAB.WEAP 

READ  11*EXPL0,REFN0,NR,RADM*R,H 

READ  12  *RHOA*CD»B*C 

READ  13*  ( (NFRAG! I  * J ) *  1  =  1  * !NT ) > J=1 *NZONF ) 

READ  14,  ( (W ( I » J ) * I  =  1  * INT ) * J=1 , NZONE ) 

READ  14, <VBAR( J) ,J=1 , NZONE) 

READ  14* ( AINT ( I ) , I =1  ,  INT  ) 

AMULT  = I MULT 

SX=H/SQRT ( H*H+4.*RADM*RADM ) 

CX=SQRT ( 1 . -SX*SX ) 

X=ARTNO(SX,CX) 

XP=X#57. 2957795131 

NZ=NZ0NE+1 

DO  100  M=2  *NZ 

THE  T A ( M ) =  THE  T  A ( M- 1 )  +  5 . 0  * AMUL T 
!'"  c>AD(M)=THETA  (Ml*. 0174  532925  1994 


Figure  A-l,  Program  Deck  Listing 


-  PRINT  26  . 

PRINT  27 

PRINT  28  . 

PRINT  29.  (AINT ( I  —  1 ) » A  INTI I ) * 1=8*1  NT) 

DO  302  J=1 .NZONE 

302  PRINT  30.  I CNBAR  (  I  ,  J  )  .CMBAR  (  I .  J )  .  I  =  8  .  INT  )  >VP,AR  ( J )  >VOC  J )  »W6AR  <  J  ) 
1N0 ( J ) 

PRINT  31.WSUN 
400  CONTINUE 
RETURN 
END 


Figure  A-l  (Cont'd) 


-  PR0=.489*RH0A*CD*R 

ANR=NR 
WSUM=0*0 
IST=C 

DO  300  J=1.NZCNE 
JJ=J+1 

DI V=SX/S IN < RAD ( JJ )  ) 

DIVI=SQRT(  l.-DIV*DI.V) 

AS  IN ( JJ) =ARTNQ(DI V.D IV1 ) 

IF(SX.GE.SIN(RAD( JJ) ).AND.SX.GE.SIN(RAD( J) ) )  GO  TO  140 
IFISX.LE.SIN (RAD(JJ) ) .AND.SX.LE.SIN(RAD( J) ) )  60-  TO  130 
I F ( IST.EO.l)  GO  TO  131  . 

AK ( J ) =THET  A ( J j ) -THETA ( J ) 

D£GX=AES ( XD-THETA ( J ) ) 

CAK( J)=2.0*OEGX  +<6.28318531*ABS(COS!RAD< JJ) )-CX) 

1* ( THETA( JJ  )-XD) )/( <RAD< JJ)-X)*( SIN (RADI JJ)  ),*ASIN( JJ )+SX*l .5707963 > 
2) 

I  ST  =  1 

SAVE=AK( J) 

SAVEC=CAK( J) 

GO  TO  142 
131  AK ( J ) =SAVE 
CAK ( J ) =SAVEC 
GO  TO  142 

140  AK ( J )  =  ( COS  <  RAD  (J  ) )-C0S<RAD< JJ) ) )*3ol4]  59265 

GO  TO  141  -  •  .  . 

1 3a  AK( J)=(SIN(RAD  ( J J ) )*AS 1 N ( J J ) +SIN ( RAD,  < J ) ) »AS INI J ) ) * ( RAD  (JJ)-RAD 
1  (  J  )  ) 

141  CAK( J)=6.28318531*<COS(RAD  (  J ) ) -COS ( RAD(J J ) ) ) 

142  CS=CAK( J)/AK(J)  • 

CAK(J)=6*28318531#( COS (RAD  (  J ) ) -COS { RAD t J J ) ) ) 

150  U=0.0 

WBAR ( J) =0.0 

SNBAR=0.0 

SME=OiO 

DO  200  1  =  1*1  NT 

'  CNBAR ( I <  J ) *CB* ( FRAG ( I » J ) /ANR ) 

CMBAR ( I » J ) =W ( I ♦ J ) /  FRAGI I » J ) 

WBAR ( J ) =CNBAR { I  *  J ) #CM8AR ( I » J ) +W6AR ( J ) 

I F  < CMBAR (1*0))  221  *221  *210 
2)0  U=U+1.0 

SMB=CM8AR ( I » J )+SMB 
221  I F ( CN8AR ( I  * J) )  200.220.220 
22u  SNBAR=CN9AR( I »J)+SNBAR 
?00  CONTINUE 

SM6AR=(1./U)*SMS 
ABAR=B*(SMBAR**(-C) ) 

ON EGA = PRO* ABA R 

•  VO ( J ) =VBAR ( J ) * ( (2.71821828459**OMEGA-l . ) /ONEGA) 

NO ( J ) =  SNBAR/CAK I J ) 

WSUM=WSUM+WBAR< J) 

300  CONTINUE 

PRINT  20.TAB.WEAP 

PRINT  21.NR.RADM.R.H.EXPL0.REFNC 

PRINT  22 

PRINT  23*AINT  (1)  » (  A INT-!  I  )  » AINT (  1  +  1  )  *1  =  1*6) 

PRINT  24 

DO  3C1  J=1»NZ0NE 

?01  PRINT  25.  THETA! J) »THETA(J+1 ). (CNBARt I » J ) * CMBAR ( I  * J ) *1=1*7) 


Figure  A-l  (Cont'd) 


Card  Type  7  , W  Total  weight  of  fragments 

per  mass- interval  per  zone 

Starting  with  the  first  polar  zone,  list  the  total  weight  of  the 
recovered  fragments  in  the  largest  mass  interval,  the  next  largest 
mass  interval,  etc.,  until  all  mass  intervals  have  been  accounted 
for.  Continue  listing. in  the  same  order  for  each  successive  polar 
zone  until  all  polar  zones  have  been  considered. 

Card  Type  8  VBAR  Average  fragment  velocity 

per  zone 

Starting  with  the  first  polar  zone,  list  the  average  fragment 
velocity,  and  continue  listing  for  each  successive  polar  zone 
until  all  have  been  considered. 

Card  Type  9  AINT  Hass  interval 

Starting  with  the  largest  mass  interval,  list  the  value  of  the 
lower  bound,  and  continue  listing  until  all  mass  intervals  have 
been  considered.  SLince  all  polar  zones  have  the  same  mass  inter¬ 
vals,  they  are  listed  only  once. 

INPUT  AND  OUTPUT 

A  sample  of  the  input  and  output  from  the  program  is  presented 
in  Figures  A-3  and  A-4.  This  table  is  based  on  data  from  static 
detonation  tests  of  four  8'/55  Mk  25  Mod  1  (HC)  projectile  warheads. 
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Figure  A-3.  Sanpla  Input. 
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Figure  A-3  (Cent'd) 
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Figure  A-4.  Sample  Output 
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Figure  A-4  (Coat'd) 


